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ABSTRACT 

£> ' We have examined the Mg-cr relation for early-type galaxies in the EFAR sample and 

0^ , its dependence on cluster properties. A comprehensive maximum likelihood treatment 

of the sample selection and measurement errors gives fits to the global Mg-cr relation of 
Mg6' = 0.131 log cr- 0.131 and Mg 2 = 0.257 log a - 0.305. The slope of these relations 
is 25% steeper than that obtained by most other authors due to the reduced bias of 
our fitting method. The intrinsic scatter in the global Mg-cr relation is estimated to 
be 0.016 mag in Mg6' and 0.023 mag in Mg2. The Mg-cr relation for cD galaxies has a 
higher zeropoint than for E and SO galaxies, implying that cDs are older and/or more 
metal-rich than other early-type galaxies with the same velocity dispersion. 

We investigate the variation in the zeropoint of the Mg-cr relation between clusters. 
Q ■ We find it is consistent with the number of galaxies observed per cluster and the 

|— | | intrinsic scatter between galaxies in the global Mg-cr relation. We find no significant 

c/) . correlation between the Mg-cr zeropoint and the cluster velocity dispersion, X-ray 

luminosity or X-ray temperature over a wide range in cluster mass. These results 
provide constraints for models of the formation of elliptical galaxies. However the Mg-cr 
relation on its own does not place strong limits on systematic errors in Fundamental 
' Plane distance estimates due to stellar population differences between clusters. 

?H \ We compare the intrinsic scatter in the Mg-cr and Fundamental Plane (FP) re- 

lations with stellar population models in order to constrain the dispersion in ages, 
metallicities and M/L ratios for early- type galaxies at fixed velocity dispersion. We 
find that variations in age alone or metallicity alone cannot explain the measured 
intrinsic scatter in both Mg-cr and the FP. We derive the joint constraints on the 
dispersion in age and metallicity implied by the scatter in the Mg-cr and FP relations 
for a simple Gaussian model. We find upper limits on the dispersions in age and metal- 
licity at fixed velocity dispersion of 32% in St/t and 38% in 8Z/Z if the variations in 
age and metallicity are uncorrelated; only strongly anti-correlated variations lead to 
significantly higher upper limits. The joint distribution of residuals from the Mg-cr and 
FP relations is only marginally consistent with a model having no correlation between 
age and metallicity, and is better-matched by a model in which age and metallicity 
variations are moderately anti-correlated (St/tm40%, 8Z/Zm50% and p«— 0.5), with 
younger galaxies being more metal-rich. 

Key words: galaxies: distances and redshifts — galaxies: elliptical and lenticular, 
cD — galaxies: stellar content — galaxies: formation — galaxies: evolution 
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1 INTRODUCTION 

The primary aim of the EFAR project (Wegner et al. 1996; 
Paper 1) is to use the tight correlations between the global 
properties of early-type galaxies embodied in the Funda- 
mental Plane (FP: Djorgovski & Davis 1987, Dressier et al. 
1987) to measure relative distances to clusters of galaxies in 
order to investigate peculiar motions and the mass distri- 
bution on large scales. However these global relations also 
constrain the dynamical properties and evolutionary histo- 
ries of early-type galaxies. For example, Renzini & Ciotti 
(1993) show that the tilt of the FP implies a range in mass- 
to- light ratio M/L among ellipticals of less than a factor of 
three, while the low scatter about the FP implies a scatter 
in M/L at any location in the plane of less than 12%. Simi- 
lar reasoning has been used to constrain the star formation 
history of cluster ellipticals using the colour-magnitude rela- 
tion (Bower et al. 1992, Kodama & Arimoto 1997). Recently 
the FP, Mg-cr and colour-magnitude relations have been 
followed out to higher redshifts and used to show that the 
early- type galaxies seen at z~l differ from present-day early- 
type galaxies in a manner consistent with passive evolution- 
ary effects (van Dokkum & Franx 1996, Ziegler & Bender 
1997, Kelson et al. 1997, Ellis et al. 1997, Kodama & Ari- 
moto 1997, Stanford et al. 1998, Kodama et al. 1998, Bender 
et al. 1998, van Dokkum et al. 1998). 

In this paper we consider the relation between the cen- 
tral velocity dispersion a and the strength of the magne- 
sium lines at a rest wavelength of 5174A for the early-type 
galaxies in the EFAR sample. This Mg-cr relation connects 
the dynamical properties of galaxy cores with their stel- 
lar populations. The remarkably small scatter about this 
relation (Burstein et al. 1988, Guzman et al. 1992, Bender 
et al. 1993, j0rgensen et al. 1996, Bender et al. 1998), and 
its distance- independent nature, make it a potentially useful 
constraint on models of the star formation history of early- 
type galaxies and a test for environmental variations in the 
FP (Burstein et al. 1988, Bender et al. 1996). 

There are, however, some problems with using the Mg-cr 
relation for probing galaxy formation. Two of these prob- 
lems are apparent from the stellar population models (e.g. 
Worthey 1994, Vazdekis et al. 1996): (i) both age and metal- 
licity contribute to the Mg linestrengths in comparable de- 
gree, so that a spread in linestrengths could be due to either 
a range of ages or a range of metallicities or some combina- 
tion; (ii) the Mg linestrengths are not particularly sensitive 
indicators — at fixed metallicity a difference in age of a fac- 
tor of ten only results in a change of 0.05-0.1 mag, while at 
fixed age a difference of 1 dex in metallicity gives a change of 
0.1-0.2 mag. Thus Mg linestrength measurements must be 
accurate in order to yield useful constraints on the ages and 
metallicities of stellar populations, and the Mg-cr relation 
on its own can only supply constraints on combinations of 
age and metallicity and not one or the other separately. 

Recently Trager (1997) has suggested that the tight- 
ness of the Mg-cr relation may be the result of a 'conspir- 
acy', in that there appears to be an anti-correlation between 
the ages and metallicities of the stellar populations in early- 
type galaxies at fixed mass which acts to reduce the scatter 
in the Mg linestrengths. Trager takes the accurate H/3, Mg 
and Fe linestrengths from Gonzalez (1993) and applies the 
stellar population models of Worthey (1994) to derive ages 



and abundances from line indices with different dependences 
on age and metallicity. He finds that at fixed velocity dis- 
persion the ages and abundances lie in a plane of almost 
constant Mg linestrength, leading him to predict little scat- 
ter in the Mg-cr relation even for large differences in age or 
metallicity — a factor of ten in age (from 1.5 Gyr to 15 Gyr) 
gives a spread in Mg2 of only 0.01-0.02 mag. This conclusion 
depends on the appropriateness of the single stellar popu- 
lation models and requires confirmation from further high- 
precision linestrength measurements. It can also be tested 
using the high-redshift samples now becoming available. 

In a similar vein, a number of authors (Ferreras et al. 
1998, Shioya & Bekki 1998, Bower et al. 1998) have recently 
re-examined whether the apparent passive evolution of the 
colour-magnitude relation out to z~l really implies a high 
redshift for the bulk of the star-formation in elliptical galax- 
ies. They conclude that in fact such evolution can be con- 
sistent with a rather broad range of ages and metallicities if 
the galaxies assembling more recently are on average more 
metal-rich than older galaxies of similar luminosity. 

As well as studies focussing on the evolution of the 
galaxy population, there have also been investigations of 
possible variations with local environment. Guzman et al. 
(1992) have suggested that there are systematic variations 
in the Mg-cr relation which affect estimates of relative dis- 
tances based on the FP. They report a significant offset 
in the zeropoint of the Mg-cr relation between galaxies 
in the core of the Coma cluster and galaxies in the clus- 
ter halo. J0rgensen and co-workers (1996, 1997) examine 
a sample of 11 clusters and find a weak correlation be- 
tween Mg linestrength and local density within the clus- 
ter which is consistent with this result. Similar offsets are 
claimed between field and cluster ellipticals by de Car- 
valho & Djorgovski (1992) and j0rgensen (1997), although 
Burstein et al. (1990) find no evidence of environmental ef- 
fects. Such systematic differences could result from different 
star-formation histories in different density environments, 
producing variations in the mass-to-light ratio of the stellar 
population. FP distance measurements would then be sub- 
ject to environment-dependent systematic errors leading to 
spurious peculiar motions. Where data for field and cluster 
ellipticals come from different sources, however, the possi- 
bility also exists that any zeropoint differences are due to 
uncertainties in the relative calibrations rather than intrin- 
sic environmental differences. 

The Mg-cr relation has thus become an important di- 
agnostic for determinations of both the star formation his- 
tory and the peculiar motions of elliptical galaxies. Here 
we examine the Mg-cr relation in the EFAR sample, which 
includes more than 500 early-type galaxies drawn from 84 
clusters spanning a wide range of environments. In §2 we 
summarise the relevant properties of the sample and the 
techniques used to determine the Mg6 and Mg2 linestrength 
indices, the central velocity dispersions cr, and the errors in 
these quantities. We present the Mg-cr relation in §3 and 
investigate how it varies from cluster to cluster within our 
sample, and with cluster velocity dispersion, X-ray luminos- 
ity and X-ray temperature. In §4 we compare our results 
with the predictions of stellar population models in order 
to derive constraints on the ages, metallicities and mass-to- 
light ratios of early-type galaxies in clusters. In particular, 
we consider the constraints on the dispersion in the ages and 
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metallicities from the intrinsic scatter in the Mg-er relation 
on its own, and in combination with the intrinsic scatter in 
the FP. Our conclusions are given in §5. 



2 THE DATA 

Here we give a short description of our sample and dataset, 
with emphasis on the velocity dispersions and line indices 
used in this paper. The interested reader can find more de- 
tail on the sample selection in Paper 1 ( Wegner et al. 1996) ; 
on the measurement, calibration and error estimation proce- 
dures for the spectroscopic parameters in Paper 2 (Wegner 
et al. 1998); and on the structural and morphological prop- 
erties of the galaxies in Paper 3 (Saglia et al. 1997). 

2.1 The sample 

The EFAR sample of galaxies comprises 736 mostly early- 
type galaxies in 84 clusters. These clusters span a range of 
richnesses and lie in two regions toward Hercules-Corona 
Borealis and Perseus-Pisces-Cetus at distances of between 
6000 km s -1 and 15000 km s -1 . In addition to this program 
sample we have also observed 52 well-known galaxies in 
Coma, Virgo and the field in order to provide a calibrat- 
ing link to previous studies. 

The EFAR galaxies are listed in Table 2 of Paper 1, 
and comprise an approximately diameter-limited sample 
of galaxies larger than about 20 arcsec with the visual 
appearance of ellipticals. Photometric imaging (Paper 3) 
shows that 8% are cDs, 12% are pure Es and 49% are 
bulge-dominated E/SOs; thus 69% of the sample are early- 
type galaxies, with the remaining 31% being spirals or 
barred galaxies. We have obtained spectroscopy for 666 
program galaxies, measuring redshifts, velocity dispersions 
and linestrength indices (Paper 2). We have used the red- 
shifts we obtained together with literature redshifts for other 
galaxies in the clusters in order to assign program galaxies 
to physical clusters. We have used the combined redshift 
data for these physical clusters to estimate cluster mean 
redshifts and velocity dispersions. The early-type galaxies 
in our sample span a wide range in luminosity, size and 
mass: they have absolute magnitudes from Ma=— 24 to 
-18 ({Mh)=-21.6; Ho^Okms -1 Mpc -1 ), effective radii 
from 1 h^, 1 kpc to 70hJ 1 kpc ((R e }=9.1 h^, 1 kpc) and cen- 
tral velocity dispersions from less than lOOkrns -1 to over 
400 km s -1 ((cr)=220kms~ 1 ). The sample is thus dominated 
by early-type galaxies with luminosities, sizes and masses 
typical of giant ellipticals. 

2.2 The measurements 

We summarise here the procedures used in measuring the 
redshifts, velocity dispersions and Mg linestrengths; full de- 
tails are given in Paper 2. 

Redshifts and velocity dispersions were measured from 
each observed galaxy spectrum using the IRAF task f xcor. 
Linestrength indices on the Lick system were determined 
using the prescription given by Gonzalez (1993). The Mg6 
and Mg2 indices were both measured: Mg2 because it is the 
index most commonly measured in previous work, and Mgb 
because it could be measured for more objects (as it requires 
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Figure 1. A summary of the errors in log a, Mgfe' and Mg2. The 
upper panel shows the cumulative distributions of the estimated 
errors, with the median and 90th-percentile errors indicated. The 
lower panel shows the calibration against the repeat observations: 
the distribution of the ratio of rms error to estimated error for 
objects with repeat measurements is compared to the predicted 
distribution assuming the estimated errors are the true errors. 

a narrower spectral range) and is better-determined (being 
less susceptible to variations in the non-linear continuum 
shape). We find it more convenient to express the 'atomic' 
Mg& index in magnitudes like the 'molecular' Mg2 index 
rather than as an equivalent width in Angstroms, since this 
puts these two indices on similar footings. The conversion is 

Mgb' =- 2 .51og 10 (l-^) , (1) 

where AA is the index bandpass (32. 5A for Mgfe). 

Error estimates for each quantity were derived from de- 
tailed Monte Carlo simulations, calibrated by comparisons 
of the estimated errors with the results obtained from re- 
peat measurements (over 40% of our sample had at least 
two spectra taken). Two sorts of corrections were applied 
to the dispersions and linestrengths: (i) an aperture correc- 
tion, based on that of j0rgensen et al. (1995), to account 
for different effective apertures sampling different parts of 
the galaxy profile, and (ii) a run correction to remove sys- 
tematic errors between different observing setups. After ap- 
plying these corrections, individual measurements for each 
galaxy were combined using a weighting scheme based on the 
estimated errors and the overall quality of the spectrum. 

The median estimated errors in the final combined val- 
ues are Aa/a = 9.1% (i.e. A log a = 0.040 dex), AMgb' = 
0.013 mag and AMg2 = 0.015 mag. The distribution of esti- 
mated errors for each quantity is shown in the upper panel 
of Figure jj]. The lower panel of the figure shows how the 
error estimates were calibrated against the repeat observa- 
tions: the distribution of the ratio of rms error to estimated 
error for objects with repeat measurements is compared to 
the predicted distribution assuming the estimated errors are 
the true errors. The initial error estimates from the simula- 
tions have been re-scaled to give the best match (under a 
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Figure 2. The Mg— a relations for all early-type galaxies with linestrength measurements in the EFAR sample: (a) The Mg&'— a relation; 
(b) the Mg2-cr relation. In both panels E galaxies are marked as circles, E/SO galaxies as squares and cD galaxies as asterisks. Errorbars 
representing the median errors in each quantity are shown in the upper right of each panel. The solid line is the ML fit and the dashed line 
is the unweighted regression of linestrength on dispersion. The distributions of residuals in Mg about the ML fit for all objects arc shown 
in the insets as open histograms. The Gaussians described by the median residual and the robustly-estimated scatter are superimposed. 
The solid histograms are the distributions of residuals for the cD galaxies, showing their offset from the overall relation. 



K-S test) to the rms errors from the repeat measurements. A 
re-scaling by factors of 0.85 and 1.15 respectively gives good 
agreement for the errors in a and Mgfe; adding 0.005 mag 
likewise gives good agreement for the errors in Mg2. 

A comparison with the literature (Paper 2, Figure 13) 
shows that our dispersions are consistent with previous 
measurements by Davies et al. (1987), Guzman (1993), 
j0rgensen (1997), Lucey et al. (1997) and Whitmore et al. 
(1985). For the subset of galaxies in common, we compared 
our linestrengths with the definitive Lick system measure- 
ments of Trager et al. (1998) in order to derive the small ze- 
ropoint corrections required to calibrate our measurements 
to the Lick system (Paper 2, Figures 14 & 15); the overlap 
of our Mg2 measurements with those of Lucey et al. (1997) 
also shows consistency (Figure 16, Paper 2). 



3 THE Mg-cr RELATION 
3.1 The global relation 

In this section we investigate the global Mg-cr relation found 
amongst the entire sample of EFAR galaxies with early-type 
morphological classifications (cD, E or E/SO; see definitions 
in Paper 3) for which we obtained linestrength measure- 
ments. The Mg6'-cr relation is shown in Figure |^a and the 
Mg2-cr relation in Figure |^b. 

In order to fit linear relations with intrinsic scatter in 
the presence of significant measurement errors in both vari- 
ables, arbitrary censoring of the dataset and a broad sam- 
ple selection function, we have developed a comprehensive 
maximum likelihood (ML) fitting procedure (Saglia et al., in 
preparation). Excluding galaxies with dispersions less than 



100 km s" 1 or selection probabilities less than 10%, and also 
outliers with low likelihoods, the ML fits to the Mg&'-cr re- 
lation (490 galaxies) and the Mg2-<r relation (423 galaxies) 
are: 

Mg6' = (0.131 ±0.017) log a- (0.131 ±0.041) , (2) 
Mg 2 = (0.257 ±0.027) log a- (0.305 ± 0.064) . (3) 

These fits are shown in Figure [| as solid lines. The ratio 
of the slopes of these relations is consistent with the Mg2- 
Mg6' relation we obtained in Paper 2: Mg 2 «1.94Mg&'-0.05. 
Monte Carlo simulations of the dataset and fitting process, 
the results of which are displayed in Figure ^, show that 
there is no bias in the ML estimates of the slopes and zero- 
points, and provide reliable estimates of the uncertainties in 
the fit. 

The ML fits can be compared to simple regressions of 
Mg6' and Mg 2 on log a. These regressions are shown in the 
figure as dashed lines, and are: 

Mg6' = (0.104 ±0.011) log cr- (0.067 ±0.026) , (4) 
Mg 2 = (0.199 ±0.016) log a- (0.168 ±0.038) . (5) 

As expected, the simple regressions yield slopes which are 
biased low due to the presence in the data of errors in the 
abscissa as well as the ordinate, and also the intrinsic scatter 
in the relation. Slightly less-biased results are obtained by 
least squares regression minimising the orthogonal residuals 
(cf. j0rgensen et al. 1996): 

Mg6' = (0.109 ±0.012) log cr- (0.078 ±0.027) , (6) 
Mg 2 = (0.215 ± 0.017) log cr - (0.205 ± 0.041) . (7) 

These least squares fits and their uncertainties are ob- 
tained using the SLOPES regression program written by 
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Figure 3. Monte Carlo simulations of the ML fits to the Mg-cr 
relations. The top panel is for Mgf/-<r and the bottom panel for 
Mg2-o\ In each panel we show the distribution of ML fits to 
the slope and zeropoint for a set of 99 simulations of the EFAR 
dataset. The input parameters (from the fits to the EFAR data) 
and the mean and standard deviation of the fitted parameters 
(from the simulations) are given in each case. 

E.D.Feigelson and described in Isobe et al. (1990) and 
Feigelson & Babu (1992). The uncertainties are under- 
estimated because these regressions do not properly account 
for the measurement errors or the selection functions. We 
conclude that previous determinations of the slope of the 
Mg-cr relation are likely to be biased low whenever the 
dataset being fitted had significant errors in the velocity 
dispersions (as is generally the case). Hereafter we adopt 
the ML fits to the Mg-cr relation. 

The distributions of the residuals in Mg6' and Mg2 
about the ML fits are shown in the insets to Figures bk 
and |^b. In order to minimise the effects of outliers, we ro- 
bustly estimate the scatter about the Mg-o" relations as half 
the range spanned by the central 68% of the data points. 
We find an observed scatter of 0.022±0.002 mag about the 
Mg&'-cr relation and 0.031±0.003 mag about the Mg2-<r re- 
lation. Excluding outliers, the distributions of residuals are 
very well fitted by Gaussians parametrised by the median 
residual and the robustly estimated scatter. There is no 
evidence for a tail of negative residuals such as noted by 
Burstein et al. (1988) and j0rgensen et al. (1996). As the lat- 
ter authors point out, the presence of such a tail is sensitive 
to the adopted slope of the Mg-cr relation. Some giant ellip- 
ticals do, however, have intrinsically weak Mg linestrengths 
for their velocity dispersions (Schweizer et al. 1990). 

The estimates of the intrinsic scatter about the rela- 
tions that are provided by the ML fit may be exaggerated 
by outliers or by deviations of the underlying distribution of 
galaxies in the Mg-cr plane from a bivariate Gaussian. We 
therefore drop the assumption of an intrinsic bivariate Gaus- 
sian distribution in the Mg-cr plane and use Monte Carlo 



0.4 - 



§• °- 3 7 

ia Z 
o 

| 0.2 - 

Q_ 

0.1 - 



0.012 0.014 0.016 0.018 0.02 
<5Mgb' (intrinsic) [mag] 



0.4 r 

£ 0.3 r 

'n _ 

o 

| 0.2 - 

o_ : 

0.1 - 



0.018 0.02 0.022 0.024 0.026 0.028 

<5Mg 2 (intrinsic) [mag] 

Figure 4. The normalised likelihood distributions for the intrin- 
sic scatter in Mgf/ and Mg2 at fixed velocity dispersion given the 
observed scatter about the Mgfe'— a and Mg2— a relations. These 
are derived from simulations based on the observed distribution 
of dispersions and linestrengths (and their errors) , assuming only 
that there is a global linear Mg-cr relation about which there is 
Gaussian intrinsic scatter. 

simulations based on the observed distribution of dispersions 
and linestrengths and their estimated errors (accounting for 
both measurement errors and run correction errors). These 
simulations assume only that there is a global linear Mg-cr 
relation about which there is Gaussian intrinsic scatter. We 
vary this intrinsic scatter and compute the robust estimate 
of the observed scatter about the fit (the half-width of the 
central 68% of the residuals) for the simulated distributions. 
The results of these simulations are presented in Figure ^, 
which shows the normalised likelihood distributions for the 
intrinsic scatter in Mgfc' and Mg2 given the observed scat- 
ter. We find that to account for the observed scatter in the 
relations we require an intrinsic scatter of 0.016±0.001 mag 
for Mgfe'-a and 0.023±0.002 mag for Mg 2 -cr. The ratio of 
the intrinsic scatter in Mg2 to the intrinsic scatter in Mgfe' 
is slightly lower than expected from the observed Mg2-Mg&' 
relation, Mg 2 «1.94Mg6'-0.05 (see Paper 2). 

Table hi compares our fits to the Mg-cr relations ob- 
tained by other authors, and gives the observed scatter 
SMgobs and the intrinsic scatter <5Mg; nt in the relations ob- 
tained in each case. For both Mg&'-cr and Mg2-cr the slopes 
we obtain are about 25% steeper than those obtained by 
most previous authors. This is not due to a difference in our 
data, but stems from our use of the ML method rather than 
regressions. In this situation regressions are biased towards 
flatter slopes than the true relation because they ignore the 
intrinsic scatter, the presence of errors in both variables and 
the selection function of the dataset. The standard or or- 
thogonal regression fits to our data, which our simulations 
show under-estimate the slope of the relations, give results 
very similar to those obtained by other authors. 

If we divide the sample by morphological type, we find 
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Figure 5. Cluster-to-cluster offsets in the Mg-cr relation. The left panels are for Mgb'-cr and the right panels for Mg2 
panels show the median offsets from the global Mg-cr relation for clusters with three or more linestrength measurements, 
panels show the distribution of offsets compared to a Gaussian with the same mean and dispersion; the scatter is 0.012 mag 
0.019 mag in Mg2. The bottom panels show the offsets as a function of redshift. 
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Table 1 . Comparison of Mg-cr relation fits 





Slope 


Intercept 


<5Mg of)s 


5M gi „ t 


Mgb' -a relation. . . 










Ziegler & Bender (1997) 


0.106 


-0.079 






EFAR (this work) 


0.131 


-0.131 


0.022 


0.016 




±0.020 


±0.048 


±0.002 


±0.001 


M32 -a relation. . . 










Burstein et al. (1988) 


0.175 


-0.110 


0.016 


0.013 


Guzman et al. (1992) 


0.260 


-0.316 


0.016 


0.013 




±0.027 


±0.003 






Bender et al. (1993) 


0.200 


-0.166 


0.025 


0.018 


j0rgensen et al. (1996) 


0.196 


-0.155 


0.025 


0.020 




±0.016 








EFAR (this work) 


0.257 


-0.305 


0.031 


0.023 




±0.028 


±0.067 


±0.003 


±0.002 



that the cDs have a zeropoint which is 0.009 mag higher than 
that of the other early-type galaxies in Mg6', and 0.014 mag 
higher in Mg2. These differences in the zeropoints are readily 
apparent in the distributions of residuals about the global 
Mg-cr relations (see the insets to Figures ^a&b), and are 
significant at the 3<r-level. Despite these zeropoint offsets, 
including or excluding the cDs changes the scatter about 
the ML fit by less than its uncertainty, as they make up 
only 10% of the whole sample. 

We find no significant differences, however, if we com- 
pare the Mg-cr relations for the two volumes of space, the 
Hercules-Corona-Borealis and Perseus-Pisces-Cetus regions, 



from which our sample is drawn. The two regions have Mg-cr 
relations with slopes and zeropoints which are consistent 
both with each other and with the overall Mg-cr relations, 
providing a check that there are no gross systematic envi- 
ronmental differences between these two regions. 

3.2 Cluster-to-cluster variations 

We do not have enough galaxies per cluster to fit both 
the slope and the zeropoint of the Mg-cr relations on a 
cluster-by-cluster basis, even in our best-sampled clusters. 
We therefore limit ourselves to investigating the variation in 
the Mg-cr zeropoint. To this end we measured the median 
offset in Mgb' and Mg2 from the global fits given above for 
the clusters with three or more linestrength measurements 
(75 clusters for Mgb and 72 for Mg2). Note that we only used 
galaxies that are cluster members based on their redshifts 
(see Paper 2). The results are not changed significantly if we 
use all clusters, or only clusters with five or more measure- 
ments. 

The top panels of Figure ^| show these zeropoint off- 
sets as a function of cluster ID number (CID), while the 
middle panels show the distributions of the offset values. 
The robustly-estimated scatter in the zeropoint offsets is 
0.012±0.002 mag in Mgb' -a and 0.019±0.004 mag in Mg 2 -cr, 
showing that the relations are remarkably uniform among 
the aggregates of galaxies in the EFAR sample. The bottom 
panels in the figure plot the same offsets as a function of red- 
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Figure 6. Variations in the Mg-cr relation with various indicators of cluster mass. The left panels are for Mgfe'-cr and the right panels 
for Mg2— cr. The top panels show the median offsets from the global Mg-cr relation as a function of log ^ ( ; the middle panels show 
the offsets as a function of log Lx ; the bottom panels show the offsets as a function of X-ray temperature, kT. 



shift, showing that there is no dependence of the relations 
on relative distance within the sample. 

This scatter in the zeropoint offsets could purely be a 
consequence of a galaxy-to-galaxy scatter in a global Mg-cr 
relation, or it could also require a variation in the zeropoint 
of the relation from cluster to cluster. These possibilities 
were examined by extending the simulations described in 
the previous section, adding a further source of scatter to the 
Mg-cr relation in the form of an intrinsic variation between 
clusters in the zeropoint of the relation. For simplicity we 
assume that this variation also has a Gaussian distribution. 

We find that if we make the extreme assumption that 
there is cluster-to-cluster scatter but no intrinsic scatter 
between galaxies within a cluster, then zeropoint varia- 
tions between clusters with an rms of 0.009 mag in Mgfe' 
and 0.015 mag in Mg2 are required to recover the ob- 
served cluster-to-cluster scatter. However this model under- 
predicts the observed scatter about the global relation, giv- 
ing 0.017±0.001 mag for Mg6' and 0.025±0.002 mag for 
Mg2 compared to the actual values of 0.022±0.002 mag and 
0.031±0.003 mag. On the other hand, if we assume that 
there is no zeropoint variation between clusters, then the 
intrinsic scatter between galaxies required to recover the 
observed scatter in the global relation (0.016 mag in Mg6' 
and 0.023 mag in Mg2; see previous section) predicts a scat- 
ter in the cluster zeropoints of 0.012±0.001 mag in Mg6' 
and 0.016±0.002 mag in Mg2, which is consistent with the 
observed values of 0.012±0.002 mag and 0.019±0.004 mag 
within the joint errors. 



We conclude that there is no evidence for significant in- 
trinsic zeropoint variations between clusters, since sampling 
a galaxy population drawn from a single global relation with 
intrinsic scatter consistent with the observations can account 
for the zeropoint differences between our clusters. 

3.3 Variation with cluster properties 

As there is very little change in the zeropoint of the Mg-cr 
relation from cluster to cluster, it follows that there can be at 
most only a weak dependence of the zeropoint on the prop- 
erties of the clusters. Here we investigate the effect of cluster 
properties on the stellar populations as reflected in the Mg-cr 
zeropoints, considering cluster velocity dispersions, X-ray lu- 
minosities and X-ray temperatures (all indicators of cluster 
mass). The cluster dispersions come from Table 7 of Paper 2, 
using redshifts both from EFAR and from the ZCAT cata- 
logue (Huchra et al. 1992; version of 1997 May 29). X-ray 
luminosities and temperatures are available for 26 of our 84 
clusters in the homogeneous and flux-limited catalogue of 
X-ray properties of Abell clusters by Ebeling et al. (1996) 
based on ROSAT All-Sky Survey data. The X-ray luminosi- 
ties are determined to a typical precision of about 20%. In 
order to have comparable precision in the cluster velocity 
dispersions, we only use clusters with dispersions computed 
from at least 20 galaxy redshifts; this also leaves 26 clusters, 
17 of which are in common with the X-ray subsample. 

Figure ^ shows the offsets in the Mg-cr relations as func- 
tions of log Clustery l°g Lx and kT. Applying the Spearman 
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rank correlation statistic, we find that there is no significant 
correlation between the Mg-cr offsets and any of these quan- 
tities, and thus no evidence for a trend in the zeropoint of 
the Mg-cr relation with cluster mass. Weighted regressions 
give best-fit relations and their uncertainties: 

AMgfo' = (-0.002 ± 0.009) logo-rf + (0.016 ±0.050), (8) 
AMg 2 = (-0.001 ± 0.011) logc4 + (0.002 ±0.060). (9) 

If we take a complementary approach, splitting the clusters 
into two subsamples about the median value of Lx and fit- 
ting the Mg-cr relations to the galaxies of the high-Lx and 
low-Lx clusters separately, we again find no significant dif- 
ferences in the slopes or the zeropoints of the fits, which are 
compatible with the global fits obtained above. 



4 DISCUSSION 

There are at least four main questions which can be ad- 
dressed using the above results. 

(i) What are the theoretical implications of the lack of 
correlation between the mass of a cluster and the zeropoint 
of the Mg-cr relation for cluster galaxies? 

(ii) What effect do the stellar population differences im- 
plied by the observed variations in the Mg-cr relation have 
on Fundamental Plane estimates of distances and peculiar 
velocities? 

(iii) What constraint does the intrinsic scatter about 
the Mg-cr relation place on the spread in age, metallicity and 
mass-to-light ratio amongst early-type galaxies in clusters? 

(iv) What further constraints on these quantities result 
from combining the scatter about the Mg-cr relation with 
the intrinsic scatter about the Fundamental Plane? 



4.1 Mg— cr zeropoint and cluster mass 

The small scatter in the zeropoint of the Mg-cr relation from 
cluster to cluster, and in particular the lack of correlation 
between the Mg-cr zeropoint and the cluster mass, seems to 
imply that the mass over-density on Mpc scales in which 
an early-type galaxy is found has little connection with its 
stellar population and star-formation history. 

The variation of the Mg-cr relation with cluster prop- 
erties has previously been studied in a sample of 11 nearby 
clusters by j0rgensen et al. (1996) and j0rgensen (1997). 
Following Guzman et al. (1992), these authors look for 
a trend in Mg2-cr offsets with the 'local density' within 
clusters. The estimator of local density used is p c i U ster = 
log o1 luster — log R, where R is the projected distance of the 
galaxy from the cluster centre. Since R is only a lower limit 
on the galaxy's true distance from the cluster centre, this is 
a rather poor estimator of the true local density. J0rgensen 
et al. find that the residuals in Mg2 show a weak trendTj with 
local density, AMg2 oc 0.009p c i U ster- Since the residuals do 
not correlate with radius within the cluster (see Figure 3 of 
j0rgensen (1997)), but do show a significant correlation with 
cluster velocity dispersion, AMg2 oc 0.02 log o-^ luster (least- 
squares fit to the data in Figure 5 of j0rgensen (1997)), we 



* The sign of the trend in equation 6 of j0rgensen (1997) is in- 
correct; the coefficient should be +0.009 (j0rgensen, priv.comm.) 



would argue that a more straightforward interpretation of 
their results is a correlation of Mg2-cr zeropoint with total 
cluster mass rather than local density. 

A correlation of this amplitude is formally consistent 
at the 2cr level with the distribution of Mg2 offsets versus 
log o\ luster for the EFAR data (see equation 0); transforming 
J0rgensen's result via the Mg2-Mg6' relation gives a corre- 
lation which is consistent at the 1.4cr level with equation [§| 
We conclude that any correlation between the Mg-cr rela- 
tion zeropoint and the cluster mass is sufficiently weak (of 
order AMg2 oc 0.02 log Cluster or less) that it is not reliably 
established by the existing data, which are consistent with 
no correlation at all. 

Semi-analytic models for the formation of elliptical 
galaxies, which previously neglected metallicity effects (see 
Kauffmann 1996, Baugh et al. 1996), are only now begin- 
ning to incorporate chemical enrichment and successfully re- 
produce the general form of the observed colour-magnitude 
and Mg-cr relations (Kauffmann & Chariot 1998). In con- 
sequence, there are as yet no reliable predictions for the 
variation of the Mg-cr relation zeropoint with cluster mass. 
The limits given above, together with limits on the difference 
in Mg-cr zeropoints for field and cluster ellipticals (Burstein 
et al. 1990, de Carvalho & Djorgovski 1992, j0rgensen 1997), 
should provide valuable additional constraints and encour- 
age further development of chemical enrichment models 
within a hierarchical framework for galaxy and cluster for- 
mation. 



4.2 Systematic effects on FP distances 

We now consider the effects on FP distance estimates of 
systematic differences in the stellar p opul ations of early-type 
galaxies from cluster to cluster. In £3.2 we found that the 



observed cluster-to-cluster variations in the Mg-cr zeropoint 
were consistent with sampling a single global Mg-cr relation 
with intrinsic scatter between galaxies, and did not require 
intrinsic scatter between clusters. Here we turn the question 
around and ask how much intrinsic cluster-to-cluster scatter 
is allowed by the observations. 



From simulations using the model described in £3.2, in- 
corporating intrinsic scatter both between galaxies and be- 
tween clusters, we find that the maximum cluster-to-cluster 
scatter allowed within the lcr uncertainties in the scatter in 
the global Mg-cr relation and the cluster zeropoints is ap- 
proximately 0.005 mag in Mg6' and 0.010 mag in Mg2. For 
our best-fit ML Mg-cr relations and a FP given by R(X(j a I f) 
with a«1.27, this level of cluster-to-cluster scatter would 
lead to rms errors in FP distance estimates of up to 10%. 
These systematic errors, resulting from differences in the 
mean stellar populations between clusters, would apply even 
to clusters in which the FP distance errors due to stellar 
population differences between galaxies had been made neg- 
ligible by observing many galaxies in the cluster. 

We emphasise that our results here do not require any 
cluster-to-cluster scatter, but are consistent with cluster-to- 
cluster scatter corresponding to systematic distance errors 
between clusters with an rms of up to 10%. We therefore 
cannot determine from the Mg-cr relation alone whether sys- 
tematic differences in the mean stellar populations between 
clusters contribute significantly (or at all) to the errors in 
FP estimates of distances and peculiar velocities. A more 
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Figure 7. The relation between Mgfe' and log M/Lr as a func- 
tion of age and metallicity in the model of Worthey (1994). 
The solid lines are contours of constant age (1.5, 2, 3, 5, 8, 12 
and 17 Gyr), with increasing line thickness indicating increasing 
age. The dashed lines are contours of constant metallicity (—0.5, 
—0.25, 0, 0.25, 0.5), with increasing line thickness indicating in- 
creasing metallicity. The dotted grid is the linear fit to the model. 



effective way of testing for such systematic differences is 
by directly comparing each cluster's zeropoint offset from 
the global Mg-cr relation to the ratio of its FP and Hubble 
distance estimates; this approach will be investigated in a 
future paper. 



4.3 Stellar population models 

To answer the questions concerning the typical age, metal- 
licity and mass-to-light ratio of early-type galaxies which 
were raised at the start of this discussion, we need to em- 
ploy stellar population models. We use the predictions from 
the single stellar population models of Worthey (1994) and 
Vazdekis et al. (1996), noting the many caveats given by 
these authors regarding their models. To simplify our anal- 
ysis, we fit Mgb', Mg2 and log M/Lr as linear functions 
of logarithmic age (logt, with t in Gyr) and metallicity 
(log Z/Zq), for galaxies with ages greater than 4 Gyr and 
metallicities in the range —0.5 to +0.5. For the model of 
Worthey (1994; Salpeter IMF) we obtain 

Mgb' as 0.0581ogi + 0.086 log Z/Zq +0.077 (10) 
Mg 2 ~ 0.1071ogi + 0.182 log Z/Zq + 0.147 (11) 
logM/L R « 0.825 log t + 0.184 log Z/Zq -0.169 . (12) 

Figure [?] compares this fit to Worthey's model in the case 
of the predicted dependence of Mgb' and log M/Lr on age 
and metallicity. The figure shows that for ages of 5 Gyr or 
greater the fit and the model are in satisfactory agreement 
for all metallicities. 

For the model of Vazdekis et al. (1996; bimodal IMF, 



^i=1.35) we have 

Mgb' » 0.051 log t + 0.083 log Z/Zq +0.084 (13) 

Mg 2 ~ 0. 115 log t + 0.187 log Z/Zq +0.137 (14) 

log M/L R » 0.673 log t+ 0.251 log Z/Zq -0.216 , (15) 

in agreement with the fit obtained by j0rgensen (1997). 
There is good agreement between the predictions of the two 
models for the dependence of Mg2 and Mgb' on age and 
metallicity, and moderately good agreement for the depen- 
dence of M/Lr. 

Note that the same change in the Mg indices is produced 
by changes in age, A logt, and metallicity, A log Z/Zq, if 
Alogt/Alog Z/Zq w 3/2. This is the '3/2 rule' of Worthey 
(1994), which applies to many of the Lick line indices, leav- 
ing them degenerate with respect to variations in age and 
metallicity. However age and metallicity produce the same 
change in log M/Lr only if A log t/A log Z/Zq ~ 1/3 or 
1/4, so that measurements of mass-to- light ratios can in 
principle be combined with Mg linestrengths to break the 
age/ met allicity degeneracy. 

4.4 Dispersion in age, metallicity and M/L 

In the following analysis we infer the dispersion in the ages 
and metallicities of early-type galaxies by comparing the 
scatter in the Mg-cr relation with the predictions of the 
single stellar population models described in the previous 
section. This analysis uses the stellar population models to 
predict differential changes in the quantities of interest, and 
not absolute values. It is also important to remember that by 
the dispersion in age or metallicity we mean the dispersion 
in these quantities at fixed log a or, equivalently, the disper- 
sion after the overall trend with log a is accounted for. Thus 
the dispersion in age or metallicity we infer is the dispersion 
at fixed galaxy mass, not the distribution of ages and metal- 
licities as a function of galaxy mass (which is related to the 
slope of the Mg-cr relation and the distribution of galaxies 
along it). 

Single stellar populations models specified by (amongst 
other parameters) a unique age and a unique metallicity can 
only provide an approximation to real galaxies, whose stel- 
lar contents must necessarily span a range (though perhaps 
a narrow one) of ages and metallicities. Since the global Mg 
indices can be quite sensitive to the detailed metallicity dis- 
tribution (Greggio 1997), some of the scatter we observe 
may be due to galaxy-to-galaxy differences in the shape of 
the metallicity distribution rather than a dispersion in the 
mean metallicity or age. 

A further complication is presented by the over- 
abundance of Mg with respect to Fe (compared to the solar 
ratio) in the cores of early-type galaxies (e.g. Peletier 1989, 
Gorgas et al. 1990, Worthey et al. 1992). As a comparison 
of Figures ^| & (?] shows, the models discussed in the previ- 
ous section are unable to account for the highest observed 
Mg linestrengths. Tantalo et al. (1998) have produced sin- 
gle stellar population models including the effects of [Mg/Fe] 
variations and find that 

AMg 2 » O.O99A[Mg/Fe]+O.O89Alogt+O.166AlogZ/Z (16) 

Comparing this equation with those above, we see that 
the differential dependence on age and metallicity is simi- 
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lar to that predicted by Worthey (1994) and Vazdekis et al. 
(1996). However, any intrinsic scatter in the [Mg/Fe]-er re- 
lation will contribute additionally to the intrinsic scatter in 
the Mg2-<r relation and reduce the dispersion in age and 
metallicity required to account for the observations. 

For these reasons, and also because of other potential 
sources of intrinsic scatter such as dark matter, rotation, 
anisotropy, projection effects and broken homology, the es- 
timates of the dispersion in age and metallicity derived here 
must be considered as upper limits. 

With these caveats in mind, we proceed to use the model 
fits given in the previous section to infer the dispersion in 
age or metallicity based on the observed intrinsic scatter of 
0.016 mag in Mgb'-a and 0.023 mag in Mg2-cr. For ease of 
interpretation we quote the dispersions in age and metal- 
licity as the fractional dispersions St/t = S log tj log e and 
SZ/Z = 8 log Z/ log e. In applying the models in what fol- 
lows, we adopt the mean of the coefficients for the two mod- 
els and give the dispersions in age and metallicity corre- 
sponding to the intrinsic scatter about the Mgb'-a relation. 
Using the intrinsic scatter obtained from the Mg2-cr rela- 
tion would give results that are ~30% smaller, since the 
observed ratio of the intrinsic scatters is <5Mg2/<5Mg&'«1.4, 
rather than about 2 as would be expected either from the 
observed Mg2-Mg6' relation or from the models. We use 
the scatter in Mgb' rather than Mg2 because our goal is to 
establish upper limits on the dispersions in age and metal- 
licity. The estimated errors in the intrinsic scatter lead to 
uncertainties in the dispersions of 5-10%. 

If age variations in single stellar populations are the 
only source of scatter then the dispersion in age is St/t— 67%, 
whereas if metallicity variations are the sole source then the 
dispersion in metallicity is SZ/Z=43%. Similarly, the ob- 
served difference in the Mg-cr relation zeropoint for the cD 
galaxies implies that these objects are either older or more 
metal-rich than normal E or E/S0 galaxies. If the zeropoint 
differences are interpreted as age differences, cDs are on av- 
erage 40% older than typical E or E/S0 galaxies (i.e. as old 
as the oldest early- type galaxies); if the zeropoint differences 
are interpreted as metallicity differences, cDs have metallic- 
ities on average 25% higher than typical E or E/S0 galaxies 
(i.e. as high as the most metal-rich early-type galaxies). 

We can also use the model fits to estimate the approx- 
imate change in M/Lr corresponding to a change in the 
Mg line indices. If these changes are caused by age varia- 
tions alone, then we find that AlogM/L_R ~ 7AMg2 and 
A\og M / Lr ~ 14AMg&'; if, however, they are due only to 
variations in metallicity we have AlogM/L_R w 1.2AMg2 
and A log M /L R ps 2.6 AMgfc' . Thus the change in log M/L R 
is about 5 times larger if the observed change in the Mg in- 
dices is due to age differences rather than metallicity differ- 
ences. The intrinsic scatter in the Mgfo'-cr relation implies a 
dispersion in mass-to-light ratio of 50% if due to age varia- 
tions, but only 10% if due to metallicity variations. 

This predicted scatter in M/L is in fact a scatter in 
luminosity or surface brightness (since that is all the models 
deal with). We can therefore readily establish the effect of 
this scatter on distances estimated using the Fundamental 
Plane (FP) if the scatter in M/L is uncorrelated with the 
galaxies' sizes and dispersions, as indeed is the case for the 
EFAR sample (at least for galaxies with <j>100 kms" ). For 
a FP given by Rota 01 ! 13 , with R the effective radius and I the 



mean surface brightness within this radius, if the scatter in 
M/L is simply a scatter in I we have A log R = /3A log M/L. 
Most determinations of the FP, including our own, yield 
f3 ~ —0.8 (e.g. Dressier et al. 1987, j0rgensen et al. 1996, 
Saglia et al. 1998). 

Combining this relation with the dependence of M/ L on 
the Mg line indices obtained above, we find that the scat- 
ter in the Mg-cr relation corresponds to an intrinsic scatter 
in relative distances estimated from the FP of 40% if due 
to age variations, or 8% if due to metallicity variations. As 
the intrinsic scatter in the FP is found to be in the range 
10-20% (Djorgovski & Davis 1987, j0rgensen et al. 1993, 
j0rgensen et al. 1996), one cannot explain both the scatter 
in the Mg-cr relation and the scatter in the FP as the re- 
sult of age variations alone or metallicity variations alone 
(unless the single stellar population models are incorrect 
or there are significant galaxy-to-galaxy differences in the 
metallicity distributions) . Suitable combinations of age vari- 
ations and metallicity variations can, however, account for 
the measured intrinsic scatter in both the Mg-cr and FP 
relations. 



4.5 Combined Mg— a and FP constraints 

As a simple model, we assume that the scatter in the FP 
and the Mg-cr relations (at fixed log cr) is entirely due to 
variations in age and metallicity (at fixed galaxy mass). 
These variations are further assumed to have Gaussian dis- 
tributions in logt and logZ/Z© with dispersions <51ogt = 
St/t loge and SlogZ = 8Z/Z loge and correlation coeffi- 
cient p (— l<p<l). While a Gaussian distribution of metal- 
licities at fixed galaxy mass is a reasonable initial hypothesis 
for describing variations in the chemical enrichment process, 
the single-peaked shape of the assumed lognormal distribu- 
tion for the mean ages may not realistically represent the 
star- formation history (even for galaxies of the same mass) . 
The dispersion in age inferred under this model should there- 
fore be considered only as a general indication of the time- 
span over which early-type galaxies of fixed mass formed the 
bulk of their stellar population. 

Writing the scatter in Mg linestrengths and FP residu- 
als as e>M g and 8-pv and the dispersion in logt and log Z/Zq 
as at and oz, this simple model relates the scatter in the 
observed quantities to the dispersion in age and metallicity 
by: 



5p P 



a 2 a 2 + 2pa t az<J t az + a 2 z a% 
b t a\ + 2pb t bzo t oz + b 2 z a% . 



(17) 
(18) 



Here a t and az are the coefficients of log t and log Z/Zq for 
Mg, and bt and bz the coefficients for log R = —0.8 log M/L, 
derived from the mean of the linear fits to the two stellar 



population models given in £4.3. 

Figure |§| shows the constraints on the variations in age 
and metallicity (assumed for now to be uncorrelated) which 
are imposed by the measured intrinsic scatter in the Mg-cr 
relations and the intrinsic dispersion in log M/ Lr inferred 
from the intrinsic scatter in the FP. The intrinsic scatter we 
find about the Mgb'-a and Mg2-cr relations is then consis- 
tent with dispersions in age and metallicity on an elliptical 
locus defined by equation (with p=0) in the St/t-SZ/Z 
plane. The different loci for Mgfo' and Mg2 (the solid lines 
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Figure 8. Constraints on the dispersion in age, St/t, and metal- 
licity, SZ/Z, from the intrinsic scatter in the Mgb'-cr relation 
(<5Mg6'=0.016 mag; thick solid line) and the Mg2-cr relation 
(<5Mg2=0.023 mag; thin solid line), and from an intrinsic scatter 
in distance about the FP of 10%, 15% and 20% (dashed lines). 



in Figure |8|) result from the difference between the observed 
ratio of the scatter in Mg2 to that in Mg6' and the predicted 
ratio from the model, and give some indication of uncertain- 
ties both in the intrinsic scatter about the Mg-u relations 
and in the model predictions. A second constraint is simi- 
larly obtained from the intrinsic scatter in distance (i.e. in 
log R) about the FP using equation [l^ (again with p—0) . 
The dashed lines in Figure ^ correspond to intrinsic scatter 
about the FP of 10%, 15% and 20%. 

The important point to note about the figure is that, as 
mentioned in §1.3, the dependences of the Mg linestrengths 



and mass-to-light ratio on age and metallicity are quite dif- 
ferent, so that (if variations in age and metallicity are un- 
correlated) the two sets of constraints are nearly orthogonal. 
Thus the region of the 8t/t—SZ/Z plane that is consistent 
with the scatter in both the Mg-<r relation and the FP is 
quite limited. If we use the intrinsic scatter in the Mgb'-<r 
relation and assume a 20% intrinsic scatter in log R about 
the FP (at the upper end of the quoted range — see, e.g., 
Djorgovski & Davis (1987) or j0rgensen et al. (1996)), we 
obtain approximate upper limits on the dispersions in age 
and metallicity of St/t— 32% and SZ/Z=38%. If, however, 
we use the intrinsic scatter in the Mg2-<r relation and adopt 




Figure 9. Constraints on the dispersion in age, St/t, and metallic- 
ity, SZ/Z, as a function of the correlation coefficient p. For values 
of p between —1 and 1 (in steps of 0.2) the constraints imposed by 
an intrinsic scatter in the Mg6'-ff relation of <5Mg6'=0.016 mag 
are shown by solid lines, while the constraints from an intrinsic 
scatter in distance about the FP of 20% are shown by dashed lines. 
The joint constraint solutions for each value of p are marked by 
circles. The uncorrelated case (p=0) is indicated by a filled circle 
and thicker lines. 



an intrinsic FP scatter of 10% (as obtained for Coma by 
j0rgensen et al. 1993), then we obtain approximate lower 
limits of St/t=15% and 5Z/Z=27%. 

Similar arguments allow us to evaluate the relative 
contributions of the dispersions in age and metallicity to 
the errors in distance estimates derived from the FP. For 
the fiducial case (oFP=20%, oMgfe'=0.016 mag and p=0), 
where St/t— 32% and 8Z/Z—3S%, the mean stellar popula- 
tion model implies that the dispersion in age gives an intrin- 
sic FP scatter of 19% while the dispersion in metallicity gives 
7%. In fact for most of the plausible range of dispersions in 
age and metallicity shown in Figure ^, it is the dispersion 
in age which dominates the intrinsic scatter about the FP. 
Only for the lowest plausible age dispersion and the highest 
plausible metallicity dispersion (St /t— 11% and SZ/Z=43%, 
corresponding to <5FP=10% and (5Mg6'=0.016 mag) does the 
contribution to the FP scatter from the dispersion in metal- 
licity achieve equality with the contribution from the dis- 
persion in age. 

The constraints on the dispersions change if there is a 
significant correlation (or anti-correlation) between the vari- 
ations in age and metallicity. Figure |^ shows how the con- 
straints corresponding to the upper limits <5FP=20% and 
SMgb' =0.016 mag (corresponding to the thick lines in Fig- 
ure ^) are modified as the correlation coefficient p varies 
over its full range from —1 to 1. Note that for p = ±1 we 
have St/t oc ^f3/2 SZ/Z . The main point to extract from 
this figure is that if the variations in age and metallicity 
have a correlation coefficient in the range — 0.5<p<l, then 
the dispersions in age and metallicity vary by only ±6% 
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Figure 10. The joint distribution of residuals from the Mgfe'-cr and FP relations: (a) the residuals for the EFAR data set; (b) a simulation 
with St/t=32% and 5Z / Z=3S% (the values derived from the intrinsic scatter in the Mgb'-tr relation and a FP scatter of 20% assuming 
the variations in age and metallicity are uncorrelated); (c) a simulation with <5t/t=67% and SZ/Z=0 (from the Mgfe'-cr intrinsic scatter 
assuming a dispersion in age only); (d) a simulation with <5t/t=0 and SZ/Z=A3% (from the Mg6'-cr intrinsic scatter assuming a dispersion 
in metallicity only). The bottom panel shows simulations consistent with the intrinsic scatter in the Mgb'-cr relation and a FP scatter of 
20%, and with various assumed correlations between age and metallicity: (e) a simulation with 6t/t=26%, SZ/Z=27% and p=l; (f) a 
simulation with <5t/i=28%, SZ/Z=31% and p=0.5; (g) a simulation with 5t/t=38%, 5Z/Z=50% and p=— 0.5; (h) a simulation with 
<54/t=57%, 5Z/Z=80% and p=— 1. The dotted lines are the expected correlations for a dispersion in age alone (A log Z=0) or metallicity 
alone (A log i=0). The dashed line is the correlation expected if the distribution is dominated by the errors in log a. 



and ±12% respectively about the values inferred in the un- 
correlated case. Only if the age and metallicity variations 
are strongly anti-correlated (pm—1; i.e. younger galaxies are 
more metal-rich) do we obtain significantly different solu- 
tions, with a broader allowed range in both age and metal- 
licity (St/t as large as 57% and 8Z/Z as large as 80%). This 
conclusion is complementary to that reached by Ferreras 
et al. (1998), who find that the apparently passive evolution 
of the colour-magnitude relation observed in high-redshift 
clusters does not necessarily imply a common epoch of ma- 
jor star-formation if younger galaxies are on average more 
metal-rich. 

We can test the degree of correlation between the vari- 
ations in age and metallicity by examining the joint distri- 
bution of residuals about the Mg&'-cr and FP relations. This 
distribution is shown for the EFAR data set in Figure p"o|a. 
There is no evidence for a correlation between the residuals 
in this figure; the Spearman rank correlation coefficient be- 
tween the residuals is 0.084, and is not significant at the 2a 
level. 



In order to investigate the expected distribution of 
residuals in the presence of the estimated measurement er- 
rors, we have performed Monte Carlo simulations of the 
EFAR data using the models for the dispersion in age and 
metallicity discussed above. Figure |To] b shows a simulation 
with <5f/t=32% and 5Z/Z=38%; these are the values de- 
rived from the intrinsic scatter in the Mg&'-cr relation and 
a FP scatter of 20% when there is no correlation between 
age and metallicity. The simulated distribution resembles 
the observed distribution, although there is a weak but sig- 
nificant anti-correlation between the residuals (due to the 
dominance of the age variations in the FP residuals) which 
is not apparent in the EFAR data. Over 100 such simula- 
tions, a two-dimensional K-S test (Press et al. 1992) gives 
a median probability of 0.3% that this distribution and the 
observed distribution are the same. 

Figures [ic|cfcd show simulated distributions for the 
cases where the intrinsic scatter in Mg&'-cr is due to age 
alone or metallicity alone. Neither case is consistent with 
the observed distribution, supporting the claim that neither 
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age nor metallicity can be solely responsible for the scat- 
ter in both the Mg&'-cr relation and the FP. Figures |ic|e-h 
show simulated distributions for four cases where the vari- 
ations in age and metallicity are correlated (with p=+l, 
+0.5, —0.5 and —1 respectively). The perfectly correlated 
and perfectly anti-correlated cases are not consistent with 
the observed distribution. However Figure |Ic| g shows that 
a distribution with no significant correlation between the 
Mg&'-cr and FP relation residuals is produced when p— — 0.5. 
A two-dimensional K-S test gives a median probability over 
100 such simulations of 1.7% that this distribution and the 
observed distribution are the same. This relatively low prob- 
ability may reflect a problem with the model, although it 
may simply be due to sampling uncertainty (the probabili- 
ties under this test vary between simulations with an rms of 
a factor of 6) or non-Gaussian outliers in the EFAR resid- 
uals. The point to be emphasised is that a model with a 
moderate degree of anti-correlation between age and metal- 
licity appears to give significantly better agreement with the 
observed distribution than a model in which age and metal- 
licity are uncorrelated. 



5 CONCLUSIONS 

We have examined the Mg-cr relation for early-type galaxies 
in the EFAR sample. We fit global Mg&'-cr and Mg2-cr rela- 
tions (equations g and ^) that have slopes about 25% steeper 
than those obtained by most previous authors. This differ- 
ence results not from the data itself but from an improved 
fitting procedure: we apply a comprehensive maximum like- 
lihood approach which correctly accounts for the biases in- 
troduced by both the sample selection function and the sig- 
nificant errors in both Mg and a. The observed scatter about 
the Mg-cr relations is 0.022 mag in Mg&' and 0.031 mag in 
Mg2; the intrinsic scatter in the global relations, estimated 
from Monte Carlo simulations, is 0.016 mag in Mg&' and 
0.023 mag in Mg 2 . 

With too few galaxies per cluster to reliably determine 
the full relation for each cluster separately, we fix the slopes 
of the relations at their global values in order to investi- 
gate the variation in the zeropoint from cluster to cluster. 
We find that the zeropoint has an observed scatter between 
clusters of 0.012 mag in Mgfo' and 0.019 mag in Mg2, and 
that this observed scatter is consistent with the small num- 
ber of galaxies sampled in each cluster being drawn from a 
single global relation with intrinsic scatter between galax- 
ies as given above — i.e. the observations do not require any 
scatter in the Mg-cr zeropoint between clusters. The allowed 
range for the intrinsic scatter between clusters corresponds 
to cluster-to-cluster systematic errors in Fundamental Plane 
distances and peculiar velocities with an rms anywhere in 
the range 0-10%. We therefore cannot determine from the 
Mg-cr relation alone whether systematic differences in the 
mean stellar populations between clusters contribute signif- 
icantly (or at all) to the errors in distances and peculiar 
velocities obtained using the Fundamental Plane. 

We have also examined the variation in the Mg-cr rela- 
tion with cluster properties. Our cluster sample ranges from 
poor clusters to clusters as rich as Coma, having velocity 
dispersions from 300kms _1 to lOOOkms -1 and X-ray lu- 
minosities spanning 0.3-8 xlO 44 ergs -1 . We do not detect a 



significant correlation of Mg-cr zeropoint with cluster veloc- 
ity dispersion, X-ray luminosity or X-ray temperature, nor 
is there any significant difference in the Mg-cr relations ob- 
tained by fitting the galaxies in the high-Lx clusters and 
low-Lx clusters separately. The predominant factor in the 
production of Mg in these early-type galaxies (and presum- 
ably other a-elements and perhaps their metallicity and 
star-formation history in general) is thus galaxy mass and 
not cluster mass. These observations place constraints on 
semi-analytic models for the formation of elliptical galaxies, 
which are now beginning to incorporate chemical enrichment 
and should soon be able to make reliable predictions for the 
variation of the Mg-cr relation with cluster mass. 

We apply the single stellar population models of 
Worthey (1994) and Vazdekis et al. (1996) to place upper 
limits on the global dispersion in the ages, metallicities and 
M/L ratios of early-type galaxies of given mass using the 
intrinsic scatter in the global Mg-cr relation. We infer an 
upper limit on the dispersion in M/Lr of 50% if the scatter 
in Mg-cr is due to age differences alone, or 10% if it is due to 
metallicity differences alone. These correspond to upper lim- 
its on the dispersion in relative galaxy distances estimated 
from the Fundamental Plane (FP) of 40% (age alone) or 8% 
(metallicity alone). Since the intrinsic scatter in the FP is 
found to be 10-20%, one cannot (within the context of the 
single stellar population models) explain both the scatter in 
the Mg-cr relation and the scatter in the FP as the result of 
age variations alone or metallicity variations alone. 

We therefore determine the joint range of dispersions 
in age and metallicity which are consistent with the mea- 
sured intrinsic scatter in both the Mg-cr and FP relations. 
For a simple model in which the galaxies have independent 
Gaussian distributions in log t and log Z/Zq , we find upper 
limits of St/t— 32% and 8Z/Z—3S% at fixed galaxy mass. 
If the variations in age and metallicity are not independent, 
but have correlation coefficient p, we find that so long as p is 
in the range —0.5 to 1 these limits on the dispersions in age 
and metallicity change by only ±6% and ±12% respectively. 
Only if the age and metallicity variations are strongly anti- 
correlated (p~— 1) do we obtain significantly higher upper 
limits, with St/t as large as 57% and 8Z/Z as large as 80%. 
The distribution of the residuals from the Mg-cr and FP 
relations is only marginally consistent with a model having 
no correlation between age and metallicity, and is better- 
matched by a model in which age and metallicity variations 
are moderately anti-correlated (8t/t~4Q%, 8Z/Z~50% and 
p~— 0.5), with younger galaxies being more metal-rich. 

Stronger bounds on the dispersion in age and metallicity 
amongst early-type galaxies of given mass will require more 
precise measurements of the deviations from the Mg-cr rela- 
tion and the Fundamental Plane and also improved models 
for the dependence of the line indices and mass-to-light ratio 
on age and metallicity. Further powerful constraints can also 
be obtained by measuring the intrinsic scatter in the Mg-cr 
and FP relations at higher redshifts, since the linestrengths 
and mass-to-light ratio have different dependences on age. 
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